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Abstract 

The distribution of the zeros of the partition function in the complex tem- 
perature plane (Fisher zeros) of the two-dimensional Q-state Potts model is 
studied for non-integer Q. On L x L self-dual lattices studied (L < 8), no 
Fisher zero lies on the unit circle po = e lS in the complex p = (eP J — \)j\fQ 
plane for Q < 1, while some of the Fisher zeros lie on the unit circle for Q > 1 
and the number of such zeros increases with increasing Q. The ferromagnetic 
and antiferromagnetic properties of the Potts model are investigated using 
the distribution of the Fisher zeros. For the Potts ferromagnet we verify the 
den Nijs formula for the thermal exponent yt- For the Potts antiferromagnet 
we also verify the Baxter conjecture for the critical temperature and present 
new results for the thermal exponents in the range < Q < 3. 

I. INTRODUCTION 

The Q-state Potts model [|I] on a lattice G for integer Q is defined by the Hamiltonian 

w=-jx;*K<7 j ) ) (i) 

where J is the coupling constant, (i, j) indicates a sum over nearest-neighbor pairs, and 
(Tj = 1, ...,Q. Fortuin and Kasteleyn have shown that the partition function can be 
written as 

Z= £ Q ni - G '\e pj - l) b{G '\ (2) 
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where the summation is taken over all subgraphs G' C G, and n{G') and b(G') are, respec- 
tively, the number of clusters and occupied bonds in G' . In Eq. (2) Q need not be an integer 
and Eq. (2) defines the partition function of the Q-state Potts model for non-integer Q. 
In this paper we discuss the ferromagnetic (FM) and antiferromagnetic (AF) properties of 
the two-dimensional Potts model for non-integer Q using the partition function zeros in the 
complex temperature plane (Fisher zeros). 

We have calculated the exact partition function of the Potts model on finite lattices of 
size 3 < L < 8 by the Chen-Hu algorithm || based on the Fortuin-Kasteleyn representation. 
We have used self-dual boundary conditions [Q for which the Fisher zeros of the Potts model 
map into themselves under the duality transformation, p — > 1/p, where p = [eP 3 — l)/y/Q. 
The self-dual lattices considered in this paper are periodic in the horizontal direction and 
there is another site above the L x L square lattice, which connects to L sites on the first 
row. 

Fig. 1 shows the Fisher zeros in the complex p plane of the Potts model for several values 
of non- integer Q. For Q < 1 no zero lies on the unit circle po = e ld and as Q approaches 
one from below some zeros approach the unit circle. In Fig. la we have omitted half of the 
Fisher zeros which give no more information because they can be obtained using the duality 
transformation from the Fisher zeros in the figure. For Q > 1 zeros begin to lie on the unit 
circle po and the number of zeros on the unit circle increases with increasing Q. Fig. lb 
shows the Fisher zeros in the range 1 < Q < 2. The two circles shown are the FM and AF 
loci for the 2-state (Ising) model in the thermodynamic limit |J . For Q > 2 the zeros in the 
AF region (Re(p) < 0) become more scattered. Finally, as shown in Fig. Id, for large Q all 
the zeros lie on the unit circle po [[§]. The distribution of zeros varies continuously with Q 
with the exception of Q = 1 and Q = 2. For Q = 1 the zeros are all degenerate at p = — 1 
||, while for Q = 2 the symmetry between the FM and AFM Ising model forces the AF 
zeros to lie on the locus p = — v^2 + e %e . 



There is no rigorous proof that the FM critical point p c — 1 for Q < 4 except for Q = 2 
For 1 < Q < 4 we observe that the zero closest to the positive real axis always lies on 
the unit circle and approaches the critical point p c = 1 in the thermodynamic limit. For 
Q < 1, however, because no zero lies on the unit circle for a finite-size lattice, we need to 
calculate the FM critical point from the zero closest to the point p c = 1 (the first zero). By 
using the Bulirsch-Stoer (BST) algorithm we extrapolated our results for finite lattices to 
infinite size. The error estimates are twice the difference between the (n — 1, 1) and (n— 1, 2) 
approximants. We find that the first zero converges on the critical point p c = 1. 

From the first zero, p±, we have calculated the thermal exponent yt{L) defined as 



II. THE FERROMAGNETIC POTTS MODEL 



Vt{L) 



ln{Abs[pi(L + 1) - l]/Abs[pi(L) - 1]} 
ln[(L + 1)/L] 



(3) 



or 



yt(L) 



ln{Im[ Pl (L + l)]/Im[ Pl (L)]} 
H(L + l)/L] 



(4) 
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For Q < 1 the imaginary part of the first zero is not a monotonic function of L, and so in 
this range we used Eq. (3) to calculate y t . For Q > 1 Eq. (4) was found to give the best 
estimate. Fig. 2 shows the BST estimates of the thermal exponent which are in excellent 
agreement with the den Nijs formula y t = (3 — 3x)/(2 — x) ||, where x = (2/7r)cos _1 (v / ^/2). 
Blote et al. || have obtained results similar to Fig. 2 using heat capacity data on infinitely 
long strips. 



III. THE ANTIFERROMAGNETIC POTTS MODEL 

For AF interaction J < the physical interval is < e^ 3 < 1, which corresponds to 
—1/y/Q < p < 0. Baxter jl(| has conjectured the existence of an AF critical point for 
Q < 3 at p c = (y/A — Q — 2)/y/Q, and we expect that in the thermodynamic limit the locus 
of the Fisher zeros cuts the negative real axis between —Xj^fQ and 0. For Q = 2.9 and 
L = 8 (Fig. lc) two zeros already lie on the negative real axis, but they are outside the 
physical interval. Table I shows the real and imaginary parts of the zero p a {L) closest to 
the AF interval of the negative real axis for Q = 0.1 and 2.9 for even-size lattices, and in 
Table I the last row are the BST estimates for L — > oo. Fig. 3 shows the BST estimates of 
the critical points of the Potts antiferromagnet from the Fisher zeros for < Q < 3 which 
are in excellent agreement with the Baxter conjecture [It]] . 



From the Fisher zeros we have also calculated the thermal exponent yt(L) defined as @ 

ln{Abs[p a (L + 2) - p c ]/Abs[ Pa (L) - p c }} 
Vt{L) = ln[(L + 2)/L] • (5) 

Fig. 4 shows the thermal exponent yt{L) for < Q < 3. In Fig. 4 as Q increases the 
thermal exponent decreases and around Q = 2 it is consistent with the known exact value 
yt — 1 of the Ising model. Note that these data are calculated from lattices of size L — 4, 6, 
and 8 only, and no attempt is made to extrapolate to infinite size. 

To our knowledge these are the first such calculations of y t for non- integer Q. Of course 
for the Ising model, Q = 2, the exact result y t = 1 has been known for many years, and 
our results are consistent with this. The value of y t for Q = 3 has been the subject of 
some debate recently. Ferreira and Sokal [|IT[] expect y t = 1/2, while Wang et al. JT2J found 



yt = 0.77. Although it might appear from Fig. 4 that our results for Q = 3 agree with 
the result of Wang et al, we must emphasize that these results are for very small lattices 
and have not been extrapolated to infinite size. If we include calculations for Q = 3 using 
the microcanonical transfer matrix for L = 10, and apply the BST algorithm our results 
are inconclusive. We are currently extending our calculations to larger lattices in order to 
address this interesting question. 
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TABLES 



TABLE I. The real part Re(p a ) and the imaginary part lm(p a ) of the zero p a {L) closest to the 
antiferromagnetic interval for Q = 0.1 and 2.9. The last row is the BST extrapolation to infinite 
size. 



Q 


0.1 




2.9 




L 


Re(> a ) 


Im(p ) 


Re(p a ) 


Im(p a ) 


4 


-0.0902065 


0.00855442 


-0.501335 


0.198694 


6 


-0.0844057 


0.00443054 


-0.507772 


0.137197 


8 


-0.0823869 


0.00272837 


-0.511662 


0.107013 


oo 


-0.078(4) 


-0.001(2) 


-0.527(5) 


0.019(41) 


Baxter's conjecture 


-0.07956 





-0.5586 
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FIG. 1. Fisher zeros in the complex p plane of 8 x 8 Q-state Potts models for (a) Q = 0.1 and 
0.9, (b) 1.1 and 1.9, (c) 2.1 and 2.9, and (d) 9999.9. In (b) and (c) two unit circles are the locus 
of the Fisher zeros for the Ising model in the thermodynamic limit. 
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FIG. 3. The critical point of the Potts antiferromagnet by the BST estimates (filled circles) 
and by the Baxter conjecture (continuous curve). 
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FIG. 4. The thermal exponent yt(L) of the Potts antiferromagnet as a function of Q for 
and L = 6. Note the large finite-size effect as Q approaches 3. 
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